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ABSTRACT 

Richness is a key defining characteristic of a galaxy cluster. We measure the optical richness of 
galaxy clusters from the CNOC1 cluster redshift survey using the galaxy-cluster center correlation 
amplitude B gc . We show that the B gc values measured using photometric catalogs are consistent 
with those derived from redshift catalogs, indicating that richness can be measured reliably from 
photometric data alone, even at moderate redshifts of ~0.6. We establish the correlations between 
optical richness and other important attributes of a galaxy cluster, such as velocity dispersion, 
mass, radius, and X-ray temperature and luminosity. We find that the scaling relations of these 
quantities with richness are entirely consistent with those derived by assuming a simple mass 
density profile at 0.5 hf 0 l Mpc of p ~ r -18 . The excellent correlations between B gc and velocity 
dispersion and X-ray temperature allow one to use richness, an easily measurable quantity using 
relatively shallow optical imaging data alone, as a predictor of these quantities at moderate 
redshifts. The B gc parameter can be used to estimate the velocity dispersion of a cluster to a 
precision of approximately 15% (~ ±100 km s _1 ), and X-ray temperature to about 20%. Similar 
correlations, but with larger scatter, are also obtained between richness and the characteristic 
radius and mass of the clusters. We compare the relative merits of B gc , Tx, and Lx as predictors 
of the dynamical mass, and find that they are comparable, providing estimates at an accuracy of 
~ 30%. We also perform similar analyses of correlations between richness and velocity dispersion, 
Tx, and Lx with a sample of low-redshift Abell clusters and find consistent results, but with larger 
scatter, which may be the result of a less homogeneous database, or sample-dependent effects. 


1. Introduction 

The single most important attribute of a galaxy 
cluster is its mass. The cluster mass function and 
its evolution provide constraints to the evolution 
of large scale structure and important cosmolog¬ 
ical parameters such as and erg (e.g., Carl- 
berg et al. 1997; Oukir & Blanchard 1997; Fan et 
al. 1997; Viana & Liddle 1999; Borgani et al. 2000; 
and many others), and possibly w^, the equation 
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of state of dark matter (e.g., Levine, Shultz, & 
White 2002). The mass-to-light ratio of clusters 
provides one of the most robust dynamical de¬ 
terminations of f l m via the Oort (1958) method 
(e.g., Carlberg et al. 1996; Girardi et al. 2000). 
Therefore, it is not surprising that over the last 
70 years, starting with Zwicky (1937), much ef¬ 
fort has been spent measuring the mass of clus¬ 
ters using a number of techniques. These include 
dynamical methods using redshift surveys (e.g., 
Kent & Gunn 1982; Carlberg et al. 1996; Gi¬ 
rardi et al. 1998, Rines et al. 2000), measur¬ 
ing the temperature of hot intracluster gas us¬ 
ing X-ray observations (e.g., Allen 1998; Lewis et 
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al. 1999), weak gravitational lensing (e.g., Hoek- 
stra et al. 1998; Clowe et al. 2000), and observa¬ 
tions of the Sunyaev-Zeldovich effect (e.g., Grego 
et al. 2001). These methods in general are ex¬ 
pensive in terms of the observational resources re¬ 
quired, especially at higher redshifts. 

There have been a number of comparisons made 
between the different methods for determining the 
mass of clusters. These include comparisons be¬ 
tween the X-ray and strong lensing methods in 
cores of clusters (e.g., Allen 1998; Wu 2001); be¬ 
tween the X-ray and weak lensing methods (e.g., 
Smail et al. 1997); and between the dynamical 
and X-ray methods (e.g., Girardi et al. 1998). 
Of particular interest is the detailed compari¬ 
son of the first three of these methods using 
the CNOC1 (Canadian Network for Observational 
Cosmology) sample of clusters performed by Lewis 
et al. (1999). Their results demonstrated that con¬ 
sistent mass estimates outside the cluster core re¬ 
gion can be obtained by the dynamical, X-ray, and 
weak lensing methods. 

Naively, in the situation where light traces 
mass, the richness of a galaxy cluster should be 
a direct indication of its mass. Furthermore, ob- 
servationally speaking, richness should be a rela¬ 
tively inexpensive parameter to measure, requir¬ 
ing only direct images of moderate depth even for 
z ~ 1 clusters. Historically, however, richness has 
not been a well-defined, quantitatively measured 
parameter. Yee & Lopez-Cruz (1999, here after 
YL99) examined in detail the properties of the 
richness parameter B gc , first used by Longair & 
Seldner (1979) for quantifying the environment of 
radio galaxies. B gc is the amplitude of the galaxy- 
cluster center correlation function measured indi¬ 
vidually for each cluster, and basically scales as 
the net counts of galaxies normalized by the lu¬ 
minosity function and spatial distribution of the 
cluster galaxies (see Seldner & Longair and YL99). 
Yee & Lopez-Cruz (LY99), using photometric cat¬ 
alogs of a large sample of low-redshift Abell clus¬ 
ters, showed that with reasonably careful atten¬ 
tion to the choice of the cluster galaxy luminos¬ 
ity function (LF) and background corrections, the 
B gc parameter is a robust measure of the richness, 
producing systematics of only 10 to 20% arising 
from the typical systematic uncertainties in the 
LF, the photometry, and the spatial distribution 
of the cluster galaxies. 


In this paper, we investigate the correlation 
of the optical richness of clusters with important 
cluster attributes such as mass, line-of-sight veloc¬ 
ity dispersion (tri), X-ray temperature (Tx) and 
luminosity (Lx), and virial radius ( r V i r ), using pri¬ 
marily the CNOC1 Cluster Redshift Survey sam¬ 
ple (Yee, Ellingson, & Carlberg 1996), currently 
the largest set of homogeneous spectroscopic and 
photometric data of galaxy clusters at the moder¬ 
ate redshift range of 0.17 to 0.55. We show that 
a well-defined and easily measurable richness pa¬ 
rameter such as B gc produces strong correlations 
with other global properties of galaxy clusters, al¬ 
lowing its use as an estimator for quantities such as 
velocity dispersion, mass, Tx, and Tx- We demon¬ 
strate that B gc and the more conventional Tx and 
Tx are good predictors of the virial mass of the 
clusters at the ~30% level, with B gc and Tx being 
somewhat superior to Tx- The correlations of the 
optical richness parameter with these key proper¬ 
ties make richness an important defining parame¬ 
ter for galaxy clusters. As we enter into an era of 
wide-field, deep optical imaging capability, more 
effective and efficient methods of finding galaxy 
clusters using optical techniques are being inves¬ 
tigated and applied (e.g., Gladders & Yee 2000; 
Kim et al. 2002). Optical cluster surveys cover¬ 
ing from tens to thousands of square degrees with 
sample size as large as 10 4 clusters are becoming 
available (e.g., the Red-Sequence Cluster Survey, 
see Yee & Gladders 2001; the SDSS, see Bahcall et 
al. 2002). The ability to use optical data alone, via 
richness, to estimate important quantities such as 
mass, velocity dispersion, and X-ray temperature 
and luminosity for clusters at large redshifts, be¬ 
comes increasingly important in the characteriza¬ 
tion of cluster samples, as obtaining X-ray or com¬ 
prehensive spectroscopic data becomes impracti¬ 
cal for such a large number of clusters. 

The remainder of the paper is organized as fol¬ 
lows. In §2, we briefly describe the optical and 
X-ray data sets. Section 3 presents the measure¬ 
ments of the richness parameter B gc for the clus¬ 
ters, including a comparison of the values deter¬ 
mined using photometric and spectroscopic data. 
Section 4 examines the correlations between the 
richness of clusters with various global properties. 
We discuss the implications and applications of 
these findings in §5, and summarize our conclu¬ 
sions in §6. Throughout this paper, unless other- 


2 



wise noted, we assume a A = 0 cosmology with 
go = 0.1, and express Hq in units of /150 = 50 km 
s -1 Mpc -1 . 

2. Cluster Samples and Data 

The CNOC1 Cluster Redslrift Survey is a 
multi-object spectroscopy survey conducted at the 
Canada-France-Hawaii 3.6m telescope (CFHT) of 
galaxies in moderately rich to rich clusters with a 
redshift range of 0.17 to 0.55. The sample consists 
of 15 high X-ray luminosity Einstein Medium Sen¬ 
sitivity Survey (EMSS) clusters (Gioia et al. 1990) 
plus Abell 2390. This sample has currently the 
most comprehensive and homogeneous optical and 
X-ray database, making it an excellent sample for 
the investigation of the correlations between clus¬ 
ter richness and other global properties. 

The observational strategy and procedure and 
the data reduction techniques for the photomet¬ 
ric and redshift catalogs of cluster galaxies are de¬ 
scribed in detail in Yee et al. (1996). The complete 
catalogs contain about 2600 redshifts, of which 
about half are considered cluster members. Ex¬ 
amples of data catalogs of the survey can be found 
in Yee et al. (1998), Ellingson et al. (1998), Abra¬ 
ham et al. (1998) and others. The velocities have 
a typical accuracy of ~135 km s -1 . Photometry 
in Gunn r and g is also available for all galaxies. 
An important feature of the CNOC1 redshift sur¬ 
vey is the estimate of statistical weights for the 
redshift catalog objects. The redshift survey used 
a sparse sampling strategy of about 1 in 2. The 
establishment of statistical weights of each galaxy 
as a function of magnitude, position, and color, al¬ 
lows the use of the redshift sample as a complete 
sample. Velocity dispersions and dynamically de¬ 
termined masses for the clusters are presented in 
Carlberg et al. (1996). The cluster MS0906+11, 
shown to be a binary in velocity space, (making 
proper dynamical analysis difficult), is excluded 
from the dynamical sample, leaving a total of 15 
clusters (Table 1). 

X-ray temperatures for the CNOC1 clusters 
were taken from the literature, primarily from 
Mushotzky & Scharf (1997) and Henry (2000). 
Sources for other temperature measurements are 
listed in Lewis et al. (1999). These tempera¬ 
tures are predominantly emission-weighted mea¬ 
sures from the Advanced Satellite for Cosmology 


and Astrophysics (ASCA), without correction for 
non-isothermality in the cluster cores. 

X-ray luminosities for these clusters are esti¬ 
mated from those calculated by Ellis & Jones 
( 2002 ) with an additonal correction for cooling 
flows. The Ellis & Jones luminosities include an 
additional aperture correction over the original lu¬ 
minosity estimates from Gioia & Luppino (1994), 
particularly for clusters with large core radii. For 
clusters which are not listed in Ellis & Jones, 
we use a similar correction to estimate aperture- 
corrected EMSS luminosities. For this correction, 
core radii were taken from Lewis et al. (1999) from 
ROSAT High Resolution Imager (HRI) observa¬ 
tions to determine the spatial profiles of the X- 
ray emission. These estimates should be generally 
consistent with those from Ellis & Jones. The lu¬ 
minosity of Abell 2390 was estimated from data 
from Pierre et al. (1996) and includes a similar 
aperture correction. We correct these EMSS lumi¬ 
nosities to bolometric luminosities using Njj val¬ 
ues from Dickey & Lockman (1990) and temper¬ 
atures derived from ASCA and Chandra observa¬ 
tions from the literature. An X-ray temperature is 
not available for MS1231+15; we assume a typical 
value of 6 keV for this correction. 

Finally, we apply a rudimentary correction for 
the effects of core cooling flows, based on the 
HRI spatial profiles from Lewis et al. (1999). For 
this correction, the central regions of clusters with 
cooling flows were omitted from beta-model fits to 
the outer profiles and the excess flux over the fit 
was attributed to the cooling flow. These correc¬ 
tions may be a slight overestimate of the effects 
of the cooling flow, however, as the spatial profiles 
were performed in the slightly softer ROSAT en¬ 
ergy band. However, these corrections are on the 
order of 2-19 percent of the total flux, and so rep¬ 
resent a small additional uncertainty. The X-ray 
temperatures and luminosities are listed in Table 
1 . 

A significant amount of data exists for lower 
redshift Abell clusters which can be used for com¬ 
parison with the CNOC1 results. In particular, 
the LOCOS (Low-redshift Cluster Optical Survey) 
sample of Lopez-Cruz & Yee (see Lopez-Cruz 1997 
and YL99) provides a large homogeneous optical 
database for 46 Abell Clusters at z ~ 0.05. The 
photometric data, based on KPNO 0.9m CCD im¬ 
ages in I c , R c , and B, were produced using the 
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same procedure as that of the CNOC1 survey. We 
have culled from the literature data on velocity 
dispersion, Tx, and Lx for clusters in this sample 
to allow us to compare the correlations of richness 
with global properties for the moderate redshift 
CNOC1 sample with a zero redshift sample. 

3. Richness Measurements of the Clusters 

We use the parameter B gc (see Longair & Seld- 
ner 1978; Andersen & Owen 1996; YL99) as an es¬ 
timate for the richness of the clusters. The B gc pa¬ 
rameter, defined as the galaxy-cluster center corre¬ 
lation amplitude (i.e., £(r) = R gc r -7 ), essentially 
measures the net number of galaxies in a cluster 
within some fixed aperture, scaled by a luminosity 
function and a spatial distribution function. 

Yee & Lopez-Cruz (YL99) have shown using 
a large sample of low-redshift Abell clusters that 
B gc is a robust richness estimate, given that the 
photometry is accurate to ~ 0.2 mag and the 
LF of cluster galaxies is known and universal to 
within about 20% accuracy in the Schechter func¬ 
tion parameters M* and a. Furthermore, they 
also showed that B gc values computed using dif¬ 
ferent absolute magnitude limits (but at least ap¬ 
proximately a magnitude past M*) and different 
sampling radii are stable to within the uncertainty. 
Yee & Lopez-Cruz also provided a calibration of 
the B gc value to the more traditional Abell Rich¬ 
ness Class. If the evolution of the cluster galaxy 
LF is known, B gc effectively allows the uniform 
computation of cluster richness for different red- 
shifts, providing a powerful parameter for the com¬ 
parison of cluster samples spanning a significant 
redshift range. 

We measure the B gc parameter for the CNOC1 
clusters via two methods: using the photometric 
catalogs and the redshift catalogs. This allows us 
to verify the accuracy of the photometric method 
which uses a statistical background count correc¬ 
tion. 

3.1. Photometric B gc 

In the standard photometric method (see, e.g., 
YL99) we count galaxies to a fixed absolute magni¬ 
tude within a standard cluster-centric radius and 
correct the expected background counts to obtain 
the net counts. We follow the prescription of YL99 
for k-correction, and count galaxies to a k- and 


evolution-corrected M r of -20.0 over a radius of 
0.5 hi o Mpc (physical). We note that by counting 
to a relatively small radius, uncertainties arising 
from projection and stochastic background varia¬ 
tions can be kept to a minimum. For background 
correction, we use counts in the r band generated 
from images of five randomly pointed fields ob¬ 
served during the CNOC1 survey. The absolute 
magnitude limit is chosen to be about 2 magni¬ 
tudes below M *, which is a relatively slowly vary¬ 
ing part of the LF, providing the most stable result 
(see YL99). 

The position of the brightest cluster galaxy 
(BCG) is chosen as the center of the cluster in the 
computation of the B gc parameter (with the BCG 
itself not included as part of the excess counts). 
We note that in most cases the BCG can be deter¬ 
mined unambiguously without redshift member¬ 
ship information by choosing the brightest galaxy 
on the red-sequence of the early-type galaxies of 
the cluster on a color-magnitude diagram. Carl- 
berg et al. (1999) show that for the CNOC1 clus¬ 
ters, the choice of using the BCG or galaxy density 
centroid as the center of the cluster produces mini¬ 
mal effect on the determination of the profile of the 
galaxy cluster. The use of the BCG as the center 
for the computation of B gc provides the advantage 
of simple consistency that can be easily applied. 

An important input in the estimate of B gc is the 
cluster galaxy LF and its evolution over the red¬ 
shift range of the sample. Ideally, the most robust 
result is obtained when the LF (and its evolution) 
is determined directly from the cluster sample it¬ 
self. However, despite being the currently largest 
cluster galaxy redshift sample at these interme¬ 
diate redshifts, the CNOC1 sample is not suffi¬ 
ciently large at the higher redshift end to deter¬ 
mine a statistically significant LF evolution mea¬ 
surement. Thus, we use the estimated CNOC2 
overall field galaxy luminosity evolution (Lin et 
al. 1999) with M*(z) ~ M(0) + Qz, where we 
adopt Q ~ 1.4 for the cluster galaxy LF. The zero 
redshift LF is taken to be a Schechter Function 
of M*( 0) = —21.9, and a = —1, which is the LF 
used by YL99 transformed to the Gunn system. 
We note that because B gc is derived using the rel¬ 
atively bright part of the galaxy LF, differential 
effects between the field and cluster galaxy LFs 
should be small. Using the three different LFs 
listed in YL99 produces less than ~ 10% varia- 
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Table 1 

CNOC1 Clusters 


Cluster 

z 

d a 
Hgc 

± 

Abell Class b 

T x (keV) 

± 

Lxbol C 

± 

MS0016+16 

0.5465 

2100 

305 

4 

8.0 

0.3 

48.00 

1.22 

MS0302+16 

0.4245 

782 

206 

0 

4.4 

0.7 

6.60 

1.56 

MS0440+02 

0.1965 

616 

170 

0 

5.3 

0.5 

7.10 

1.50 

MS0451+02 

0.2011 

1275 

222 

2 

8.6 

0.3 

15.47 

3.30 

MS0451-03 

0.5391 

2060 

301 

4 

10.4 

0.75 

55.40 

2.64 

MS0839+29 

0.1930 

1318 

223 

2 

4.2 

0.2 

7.87 

1.49 

MS1006+12 

0.2604 

1474 

239 

2 

8.5 

1.0 

13.26 

3.10 

MS1008-12 

0.3063 

1757 

258 

3 

7.3 

1.0 

14.97 

3.40 

MS1224+20 

0.3255 

508 

176 

0 

4.3 

0.65 

5.86 

1.34 

MS1231+15 

0.2353 

904 

196 

1 

— 

— 

11.85 

2.68 

MS1358+62 

0.3290 

1373 

238 

2 

6.5 

0.3 

14.85 

0.66 

MS1455+22 

0.2568 

593 

176 

0 

5.5 

0.15 

24.28 

2.02 

MS1512+36 

0.3717 

610 

189 

0 

3.6 

0.6 

6.17 

1.30 

MS1621+26 

0.4275 

1158 

231 

1 

6.6 

0.9 

18.90 

1.85 

A2390 

0.2280 

1604 

243 

3 

8.9 

0.35 

45.43 

14.0 


a B gc in units of /ijg 1 ' 8 Mpc 1 ' 8 . 
b Abell richness class estimated from B gc . 
c Bolometric X-ray luminosity in units of 10 44 /ii~ 0 2 erg s 
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tions in the B gc values. Furthermore, adjusting Q 
by ±30% does not change the relative B gc values 
of the high and low redshift clusters significantly. 

The resultant B gc measurements are listed in 
Table 1, along with their estimated Abell rich¬ 
ness classes. The B gc values, in units of hV A ]' 
Mpc 1 ' 8 can be roughly calibrated to the Abell 
richness class scale as follows (YL99): Abell 0: 
400 < B gc < 800; Abell 1: 800 < B gc < 1200; 
and Abell 2: 1200 < B gc < 1600. 
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Fig. 1.— Photometric B gc vs spectroscopic B gc , 
showing that richness can be estimated robustly 
using photometric data alone. The dashed line 
has a slope of 1. 

3.2. Spectroscopic B gc 

The CNOC1 sample offers the unique oppor¬ 
tunity to verify the robustness of photometrically 
determined richness parameters with those deter¬ 
mined using a redshift sample. This allows us to 
test the validity and accuracy of statistical back¬ 
ground corrections. We determine the B gc param¬ 
eters for the clusters using the same basic proce¬ 
dure as that for the photometric measurements. 
Here, galaxies considered to be in the cluster, 
based on the redshift range established for each 
cluster in Carlberg et al. (1996), are summed ac¬ 
cording to their magnitude and geometric weights 
(Yee et al. 1996) to arrive at the estimated net 


counts. In general the spectroscopic B gc 's are de¬ 
termined using absolute magnitude limits brighter 
than -20 (due to the relatively bright spectro¬ 
scopic limits), and hence have larger uncertain¬ 
ties despite not having to correct for background 
counts. 

The B gc values derived from the spectroscopic 
sample are plotted versus the photometric B gc in 
Figure 1. There is a very strong correlation be¬ 
tween the two measurements with a slope close 
to the expected unity. The scatter in the rela¬ 
tion is consistent with the estimated uncertain¬ 
ties. From this we can conclude that B gc values 
determined using photometric data with statisti¬ 
cal background count corrections are reliable es¬ 
timates of cluster richness, even to redshifts of 
~ 0.6. This further demonstrates that photomet¬ 
ric B gc is a robust richness parameter. 

4. Correlations of Richness and Cluster 
Global Properties 

4.1. Expected Scaling Relations 

The richness parameter B gc scales with the 
net counts of galaxies to some fixed luminosity 
within a certain radius. By making the simpli¬ 
fying assumption that galaxy clusters have similar 
galaxy luminosity functions at the bright end (e.g., 
Lopez-Cruz 1997 found that M* for the LOCOS 
sample has a scatter of < 0.25 mag), galaxy spatial 
distributions, and mass-to-light ratios (e.g., Carl¬ 
berg et al. 1996), B gc is then expected to track the 
mass of clusters, and essentially measures the scal¬ 
ing of the spatial distribution of galaxies around 
the center of the cluster. Hence, B gc scales as the 
mass of clusters at a fixed proper radius. How¬ 
ever, this is not an interesting correlation, as mass 
measures such as M V i r , the virial mass, or Af 200 > 
(the mass interior to r 200 j where the average mass 
density is 200p c ), are connected more directly to 
observational quantities such as o\ and lx- We 
can derive the links between B gc and these prop¬ 
erties using the following scaling relations. 

Assuming a correlation function form of £(r) = 
B gc r “ 7 and that galaxies trace the total mass, we 
expect the mass density p ~ r -7 , and hence the 
mass interior to some radius r, M(r) oc f? gc r -7 . 
The mass of a cluster at some fixed fiducial ra¬ 
dius r' is related to the virial mass M v i r at ra¬ 
dius r V i r by M r i oc M V i r r^,~ 3 . In the virialized 
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region, r V i r oc u\ (e.g., see Carlberg et al. 1996) 
and M vir oc erf. Hence, we expect the simple rela¬ 
tionship: 

B gc cx <tJ, (1) 

and 

B gc cx Mjl 3 . (2) 

Note that the same relationship holds for M 2 oo- 
Since for a virialized system M V i r oc T^ 2 (e.g., 
see Horner et al. 1999), we also expect: 

Bgc oc Tl' 2 . (3) 

Finally, another important parameter for clusters 
is r 2 oo, which provides one with a physical size 
scaling of clusters of various richness. Since r 2 oo oc 
cri , we have: 

Bgc oc rcJoo- (4) 


o 

o 

CD 


O 


O 

O 

Cci 


ID 

ft 


O 

O 



1000 2000 3000 


We note that for computing B gc we have cho¬ 
sen to use 7 = 1 . 8 , primarily to allow for a di¬ 
rect comparison with the galaxy-galaxy correla¬ 
tion amplitude. Analyses of correlations between 
galaxies and clusters have shown that 7 is typi¬ 
cally steeper than that for galaxy-galaxy correla¬ 
tion functions, varying from 1.8 to 2.4 (e.g., Lilje 
& Efstathiou 1988; Moore et al. 1994; Croft, Dal¬ 
ton, & Efstathiou 1999), with the low value corre¬ 
sponding to the correlation between clusters and 
IRAS galaxies. However, for the purpose of using 
the correlation amplitude as a measure of richness, 
a shallower 7 of 1.8 is likely the more appropriate 
one, due to the small cluster-centric radius used 
for the counting of galaxies. The 7 ’s derived for 
galaxy-cluster correlations are generally measured 
over scales of 10 h^ 0 Mpc or more, compared to 
0.5 hj q 1 Mpc used for the estimate of B gc . These 
correlation function slopes essentially reflect the 
dependence of cluster profiles at large radii. Pro¬ 
files such as NFW (Navarro, Frenk, & White 1997) 
and Hernquist models (Hernquist 1990) typically 
have a form of r <-2 at large r, and asymptotically 
approach r <_ , similar to the measured correla¬ 
tion functions. At smaller radii, however, cluster 
profiles in general have flattening slopes, possibly 
leading to a core. At the typical 0.5 h^ 0 Mpc 
where B gc is measured (equivalent to r\j 0.25r 2O o), 
a 7 of 1.8 is a reasonable representation. 


Richness (B„) 

Fig. 2.— Line-of-sight velocity dispersion vs 
richness as parametrized by B gc . The outlying 
point on the left hand edge (open triangle) is 
MS 1445+22. The dashed line is the best linear fit 
with MS 1445+22 removed, while the dotted line 
is the power-law fit. 

4.2. Correlations of Richness with Global 
Optical Properties 

In this section we examine the correlation of 
the richness of the CNOC1 clusters with quanti¬ 
ties derived from optical data, such as velocity dis¬ 
persion, mass, and r 2 oo- The main motivation in 
deriving these dependences is to use B gc , an eas¬ 
ily measurable parameter using just imaging data, 
as a predictor of the other quantities. Hence, we 
fit expressions of these quantities as a function of 
Bgc- 

The B gc parameter provides an excellent cor¬ 
relation with velocity dispersion. We plot in Fig¬ 
ure 2 the correlation of velocity dispersion (from 
Carlberg, Yee, & Ellingson 1997) vs the photo¬ 
metric B gc richness parameter using linear scal¬ 
ing, and in Figure 3 using logarithmic scaling. In 
the plots, there is one cluster which stands out 
as an outlier: MS 1455+22, which has too large 
a velocity dispersion for its richness. We exclude 
MS 1455+22 from all our fits, but include it in the 
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Fig. 3.— Velocity dispersion vs B gc in logarithmic 
scale. The dashed line is the best fitting function 
using BCES. MS1455+22, omitted from the fit, is 
indicated by an open triangle. 

plots. We further discuss the case of MS 1455+22 
in §5. Using the bisector BCES estimator fitting 
routine (which allows for measurement errors in 
both variables; Akritas & Bershady 1996) in loga¬ 
rithmic space, we obtain a fit plotted as a dashed 
line in Figure 3: 

logcri = (0.55 ± 0.09) log B gc + (1.26 + 0.30), (5) 

where o\ is in units of km s _1 , and B gCl in units 
of /bjk ' 8 Mpc 1-8 . The fit is equivalent to B gc oc 
<t(’ 8=f °' , entirely consistent with the expected re¬ 
lation of eqn 1 for 7 = 1.8. We note that varying 
7 in the derivation of B gc produces values that are 
simply scaled by the integration constant J 7 (de¬ 
fined in Longair & Seldner 1979), and does not af¬ 
fect the power-law slope of the fit. Hence, the cor¬ 
relation between velocity dispersion and B gc indi¬ 
cates that the 7 used is consistent with the galaxy 
density slope near 0.5h^ 0 Mpc. 

Since the data span a factor of only ~ 3 in rich¬ 
ness, we also perform a linear fit over the data 
range to provide a convenient formula for estimat¬ 
ing the velocity dispersion from B gc : 

ai = (416 ±61)^gj + (340 + 90) km s"\ ( 6 ) 


in the range of 500 < B gc < 2500. In Figure 2, we 
plot both the linear and power-law fits between ay 
and B gc . 

In Figures 4 and 5, we plot in logarithmic scale 
the relationship between B gc and M 200 and r 2 oo> 
with the following best BCES fits: 

log M 200 = (1.64 + 0.28) \ogB gc + (10.05 + 0.89), 

(7) 

and 

logr 2 oo = (0.47+0.16) logR sc -(1.05 + 0.48), ( 8 ) 

where M 200 is in units of and r 200 i hr 

units of h^ 0 Mpc. These correlations are equiva¬ 
lent to B gc oc Af 2 QQ 1±al °, and B gc oc rf^ 0 ' 7 , in 
excellent agreement with the expected values of 
0.60 and 1.8 for 7 = 1.8. We note that the scatter 
in the r 2 oo versus B gc correlation is rather large, 
and significantly greater than the error bars. 



Fig. 4.— M 200 vs B gc . The dashed line shows the 
best fitting power law. MS 1455+22 is indicated as 
an open triangle. 

4.3. Correlations of Richness with X-ray 
Properties 

X-ray temperature and luminosity have been 
used as key measures of cluster properties. From 
simple physical models, it is expected that these 
quantities scale with properties such as mass and 






Fig. 5.— ?’ 2 oo vs B gc . The dashed line shows the 
best fitting power law. MS 1455+22 is indicated as 
an open triangle. 

velocity dispersion. Given the excellent correla¬ 
tions of the richness parameter B gc with mass and 
velocity dispersion, the B gc parameter should also 
be a useful predictor for X-ray properties. 

Figure 6 presents the correlation between 
Tx and B gc . The bisector fit for the relation is: 

log lx = (0.78 ± 0.13) log B gc - (1.59 + 0.40), (9) 

where Tx is in units of keV. Several of the EMSS 
clusters have their total fluxes affected by cooling 
flows in the cluster core (Lewis et al. 1999), which 
may affect the temperature measurement. The 
four clusters expected to be most affected by cool¬ 
ing flows (MS0839+29, MS1224+20, MS1358+62, 
and MS1445+22) are marked by separate sym¬ 
bols on Figure 6 . We repeat the fit with these 
objects removed (leaving 10 clusters), producing 
essentially an identical result. The Ty^-B gc scal¬ 
ing of Eqn (9) corresponds to B gc oc j£- 28±0 - 21 ; 
almost 2 sigma higher than the expected power- 
law exponent of 0.90 for 7 = 1.8. The expected 
relationship, however, is predicated on 07 oc T-£- 5 . 
Different investigators have obtained different de¬ 
pendences between 07 and Tx, with exponent val¬ 
ues of 0.62 + 0.04 (Girardi et al. 1998), 0.67 + 0.09 


(Wu et al. 1999), and 0.53 ± 0.04 (Horner et 
al. 1999). Figure 7 shows the 07 vs Tx relation 
for the CNOC1 sample. The best fit gives a rela¬ 
tionship of 07 a j£- 66±0 - 13 ; closer to the steeper 
of the current literature results. Assuming such 
a Tx-cr 1 relationship, we would expect a steeper 
correlation, B gc oc +x , which is well within one 
sigma of the result obtained. 


C\2 



Fig. 6 .— X-ray temperature vs B gc . The dashed 
line shows the best fitting power law using all data 
except MS 1455+22. MS 1455+22 is indicated as 
a triangle. Three other clusters with prominent 
cooling flows are indicated by imbedded solid cir¬ 
cles. These three clusters, along with MS1455+22, 
produce a significant fraction of the scatter in the 
relation. Note that the most discrepant point is 
not MS1455+22, but MS0839+29, which is the 
point well below the fitted relation. 

The relationship between the bolometric Tx and 
B gc is plotted in Figure 8 . The best BCES fit is: 

logT x = (1.84 + 0.24) log-Bg C — (4.48 + 0.75), (10) 

where Tx is in units of 10 44 /i ^" 0 2 erg s _1 . The cor¬ 
relation is equivalent to B gc oc T^ 54±0 ' 07 . We note 
that careful attention must be given to the various 
corrections required for the determination of the 
bolometric Tx, which may introduce a change of 
a factor two or more to the result (e.g., see Ellis & 
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Log T x (keV) 



Fig. 7.— Velocity dispersion vs Tx- The dashed 
line shows the best fitting power law. The dot¬ 
ted line has a slope of 0.5, the expected scaling 
between Tx and a. 


Fig. 8.— Lx vs B gc . The dashed line shows the 
best fitting power law. MS1455+22 is indicated 
by an open triangle. 


Jones 2002). A considerably larger scatter in the 
correlation is obtained using existing values from 
the literature (e.g., from the compilation of Wu et 
al. 1999; or from the EMSS measures from Gioia 
& Luppino 1994). 

5. Discussion 

5.1. B gc as an Estimator of Cluster Prop¬ 
erties 

We have found that for the CNOC1 clusters 
strong correlations exist between the optical rich¬ 
ness and other directly measured quantities such 
as velocity dispersion and Tx, especially after a 
single discrepant cluster, MS1455+22, is removed. 
Here, we examine the goodness of fit tests for the 
various relationship to obtain an indication of how 
well B gc can be used as a predictor of other global 
properties of clusters. 

The rms scatter of the CNOC1 sample veloc¬ 
ity dispersions from the best fitting power-law re¬ 
lation, excluding MS1455+22, is 105 km s -1 , or 
~12% of the mean value. The reduced % 2 of 
the fit is 0.65. We note that, because there are 
uncertainties in the data in both the x and y 


axes, x 2 is estimated using the effective variance: 
<= cry + a 2 cx 2 , where a is the slope of the 
power-law fit. By including MS 1455+22, the rms 
deviation is increased to 176 km s -1 (or 20%), 
with the reduced y 2 increased to 1.3. Similar re¬ 
sults are obtained when using the linear relation 
of Eqn 6, but with larger reduced y 2 of 1.4 and 
2.4 for with and without MS1455+22, respectively. 
These somewhat higher y 2 values are probably a 
more robust indication of the goodness of fit, as 
the effective variance used is strictly correct only 
for a linear fit. The relatively small % 2 and frac¬ 
tional scatter suggests that, with the exception of 
MS 1455+22, the a\ vs richness relation is consis¬ 
tent with little or no intrinsic scatter beyond the 
measurement uncertainties, and B gc can be used 
as a reasonably accurate estimator of velocity dis¬ 
persion. 

For the dynamically derived quantities, such 
as mass and r 2 oo> the scatter of the data points 
around the best fitting relation is significantly 
larger than that of a\. The rms scatter in M 200 
from the best fitting power law is 4.4x 10 14 ft^ o 1 MQ 
amounting to 31% of the average value of the sam¬ 
ple; while for r 2 oo> the scatter is 0.52 Mpc, or 
21% of the average value. The reduced y 2 ’s are 
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0.83 and 2.0 for the AT 2 oo and r 2 oo relation, re¬ 
spectively. 

In general the correlations between richness and 
X-ray quantities have somewhat larger x 2 ’s. For 
the 7x vs B gc relationship, we obtain for the 13 
data points (excluding MS1455+22) a scatter of 
1.36 keV relative to the best-fit power law, about 
21% of the mean Tx of the sample. However, a 
relatively large reduced y 2 of ~ 2.1 is obtained 
for the best fit, indicating that there may be a 
significant intrinsic variance in the Tx vs richness 
relationship. Or alternatively, this scatter could 
be an indication that there may be systematic un¬ 
certainties in the Tx determinations not taken into 
account by the error bars listed in the literature. 
The correction of Tx for cooling flows and depar¬ 
tures from isothermal conditions is a complicating 
issue (e.g., see Allen 1998; Lewis et al. 1999), 
and may contribute to the larger %• Removing 
the three additional clusters with the most promi¬ 
nent cooling flows, a scatter of 1.2 keV about the 
best fitting relation with a reduced % 2 of 1.2 is 
obtained, indicating that cooling flows in clusters 
may be an important source of the larger scatter 
seen in this relationship. 

The correlation between Lx and B gc has 
the largest scatter. The reduced y 2 (excluding 
MS1455+22) of the best fit is ~ 1.5, and the rms 
scatter of the data from the fit is 47%. Removing 
the three additional cooling flow prominent clus¬ 
ters improves the goodness of fit only marginally. 

5.2. B gc , Tx, and Lx as Estimators of Dy¬ 
namical Mass 

Both Tx and Tx have often been used as esti¬ 
mators or proxies for the mass of clusters in the 
study of mass functions of clusters (e.g., Oukir 
& Blanchard 1997; Viana & Liddle 1999; Henry 
2000). Hence, it is useful to examine the corre¬ 
lations between Tx and Tx with the dynamically 
derived A7 2 oo of this sample, and compare their 
relative merits as estimators of mass to that of 
Bg C - 

Figure 9 shows the excellent correlation be¬ 
tween Tx and Af 2 00, with the best power-law fit: 

log M 200 = (2.01 ± 0.33) logTx + (13.47 + 0.29), 

(ID 

where M 2 qo is in units of h 5 Q M@, and Tx in units 
of keV. We note that MS 1455+22 is a significant 


outlier also in this correlation - it has too high a 
virial mass for its Tx- For simple virialized sys¬ 
tems, we expect AT 20 o oc T^' 5 . The relation de¬ 
rived here has a higher exponent, but is marginally 
consistent at 1.5 sigma with the scaling expected 
for virialized systems. The result, however, is con¬ 
sistent with a number of studies of X-ray clusters 
showing that M oc T^ 17 (e.g., Ettori & Fabian 
1999; Nevalainen, Markevitch, & Forman 2000). 
We note that these studies use mass determined 
using X-ray data, whereas the CNOC1 relation is 
for dynamically determined mass. The similar re¬ 
sults indicate the consistency in the two mass de¬ 
termination techniques. 

We find Tx also correlates well with M 2 00 , as 
shown in Figure 10, with: 

log M 200 = (0.89 ±0.11) logTx + (14.03 + 0.14), 

( 12 ) 

where Tx is in units of 10 44 hj erg s A The 
correlation is consistent with a linear relationship 
between Tx and mass. The self-similar isother¬ 
mal model with pure bremsstrahlung radiation 
and viral equilibrium leads to the scaling relation 
of Tx oc T^ 5 , (e.g., see Arnaud & Evrarcl 1999). 
Assuming the relation from simple virialized sys¬ 
tems, M 20 0 oc T^ 5 , we would expect AT 20 o oc T^ 75 , 
which is within 1.5 sigma of the relation obtained 
here with the CNOC1 sample. If we use the 
measured M 200 oc Tx' 01±0 ' 33 , we would predict 
A7 20 o oc Tx° ±0 ' 17 , in good agreement with the 
measured correlation. In investigations of avail¬ 
able X-ray data of galaxy clusters, it has been 
found in general that Tx oc T^ 0 " 5 (e.g., Markevitch 
1998; Allen & Fabian 1998; Arnaud & Evrard 
1999; see Babul et al. 2002 for detailed modeling 
of the Tx-Tx relation). The current data set pro¬ 
duces a scaling of Tx oc T^ 42±0 ' 06 , consistent with 
the correlations found for high temperature sys¬ 
tems (e.g., Allen, Schmidt, & Fabian 2001; Babul 
et al. 2002). 

The CNOC1 sample has shown that all three 
quantities: richness, x-ray temperature, and x-ray 
luminosity, correlate well with the derived dynam¬ 
ical mass. Here, we compare the relative merits of 
these quantities as predictors of the mass of galaxy 
clusters. The B gc parameter produces a scatter in 
the predicted Af 20 o around 31%, with a reduced 
X 2 of 0.83. This compares with a scatter in A7 20 o 
of 29% and a reduced y 2 of 1.6 for the Tx-Af 20 o 
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relation, and a scatter of 35% and a reduced y 2 
of 2.8 for Lx, all computed without MS1455+22. 
Thus, all three quantities can serve as good pre¬ 
dictor of the cluster dynamical mass to the level 
of ~ 30%, with Tx and B gc being somewhat supe¬ 
rior to Lx- However, we note that the proper de¬ 
termination of the bolometric Lx is not a simple 
procedure, and in fact requires the knowledge of 
Tx and the X-ray emission spatial profile. X-ray 
luminosities derived without careful attention to 
spatial and spectral distributions produce a much 
poorer correlation with mass and other quantities 
than shown here. Hence, Lx is not as useful an 
estimator of cluster mass as Tx or B gc . 



Fig. 9.— Dynamically determined mass M 200 vs 
Tx- The dashed line shows the MS 1455+22 is in¬ 
dicated by an open triangle. 

That optical richness can predict the mass of 
clusters as accurately as X-ray properties is not 
entirely surprising. Carlberg et al. (1996) showed 
that for the CNOC1 sample L200, the luminosity 
within 7*200, is well correlated with velocity disper¬ 
sion. Girardi et al. (2002), using a large sample of 
low-redshift clusters, also found a good correlation 
between the B band luminosity within the viral 
radius and virial mass with M v oc L^ 1 ' 2 . How¬ 
ever, neither of these luminosity measurements is a 
simple quantity to derive directly from photomet¬ 
ric data alone, as it requires the knowledge of the 
virial radius or 7 * 200 - The derivation of B gc does 
not require the knowledge of the size of the clus- 



Fig. 10.— Dynamically determined mass M 200 vs 
Lx- The dashed line shows the best fitting power 
law. MS1455+22 is indicated by an open triangle. 

ter because it assumes a fixed form of the density 
function. Furthermore, such an assumption allows 
one to estimate the characteristic size of the clus¬ 
ter (e.g., 7*200)- The correlation between B gc and 
dynamical mass derived here spans only a small 
mass range of about a factor of 10. There is a 
significant amount of both observational and the¬ 
oretical evidence that, at least at lower masses, 
the assumption of a constant mass-to-light ratio is 
not a valid one (e.g., Carlberg et al. 2001; Girardi 
et al. 2002; Balrcall et al. 2000). This will affect 
the H gc -mass correlation, and thus, the scaling 
relations between B gc and 0 * 1 , M200, and 7*200 de¬ 
rived in §4 probably cannot be extrapolated into a 
mass regime much smaller than 10 14 Mq. Further 
calibrations of these relations are needed using a 
larger mass range. 

5.3. The Case of MS1455+22 

The galaxy cluster MS1455+22 is a significant 
outlier to most of the correlations discussed in §4: 
It has either too low a B gc value (by about a factor 
of three) for its velocity dispersion, too high a ve¬ 
locity dispersion (by about a factor of two) for its 
richness, or a combination of both. In the distribu¬ 
tion of mass-to-light ratios of the CNOC1 clusters, 
MS1455+22 also stands out with an anomalously 
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large (at the 2er level) value of 807 h±oo M 0 /L 0 , 
compared to the mean of 295 hioo M 0 /L 0 for the 
whole sample (Carlberg et al. 1996). In many ways 
it is not surprising to find outliers to the correla¬ 
tions. Galaxy clusters are complex systems; differ¬ 
ential evolution and events such as major mergers 
should produce intrinsic variances to these rela¬ 
tions. Hence, it is of interest to examine clusters 
that are significantly discrepant from the average 
population. 

MS 1455+22 itself appears to be a very evolved 
cluster with a luminous and large central galaxy, 
a well-formed red sequence with a relatively small 
blue fraction (0.15; Ellingson et al. 2000), and reg¬ 
ular X-ray morphology (Lewis et al. 1999). In ad¬ 
dition, it has one of the largest cooling flow mea¬ 
surements in the X-ray (Lewis et al. 1999). The 
fact that the spectroscopic and photometric 73 gc ’s 
agree within the uncertainties indicates that the 
B gc value is not affected by an unusually low back¬ 
ground count. If B gc is anomalously low compared 
to other clusters, one possible explanation is that 
MS 1455+22 is in a more advanced state of evo¬ 
lution than the others, as suggested by the large 
cooling flows, the strong red sequence, and the 
dominant central galaxy. The dominant central 
galaxy of MS1455+22 could contribute to the low 
B gc value by having cannibalized a larger number 
of its bright members, leaving a lower net cluster 
galaxy counts. The more evolved galaxy popula¬ 
tion could produce a lower average luminosity for 
the galaxies, which again may produce a smaller 
B gc value. However, the small number statistics 
of the redshift sample are not sufficient for us to 
draw definitive conclusions regarding the LF of the 
galaxies. 

The velocity dispersion for MS1455+22 was de¬ 
termined using 51 cluster galaxies (Carlberg et 
al. 1996), and hence is expected to be relatively 
robust. Nevertheless, there are indications that 
the 1196 km s -1 value could be somewhat high. 
One indication can be found in the a\ vs Tx plot 
where MS 1455+22 has a high velocity dispersion 
for its Tx, while in the Tx versus B gc correlation 
MS 1455+22 is not particularly discrepant. How¬ 
ever, a typical cooling flow correction may bring 
Tx to a better agreement with the velocity disper¬ 
sion. Second, Borgani et al. (2000) applied several 
different methods to the CNOC1 data in deter¬ 
mining the velocity dispersion, showing that the 


results, as a sample, are consistent. However, the 
different methods in Borgani et al. produced con¬ 
sistently lower <7i for MS1455+22, ranging from 
962 km s -1 to 1033 km s -1 . (We note that us¬ 
ing the different velocity dispersion derivations in 
Borgani et al. produce essentially the same correla¬ 
tions with Bgc, but with somewhat larger scatter, 
with MS1455+22 being a less discrepant point.) 

Finally, in comparing the masses determined 
using X-ray data and weak lensing, Lewis et 
al. (1999) found that while the dynamical mass 
of MS 1455+22 is consistent with the weak-lensing 
mass, it is 1.34 times that of the X-ray mass. 
Since weak lensing tends to over estimate cluster 
mass on the average (Lewis et al. 1999), this may 
suggest that the dynamical mass of MS 1455+22 
is somewhat high. 

In summary, current data do not provide a 
single clear explanation of the discrepancy of 
MS 1455+22 in an otherwise tight correlation of <J\ 
with B gc . There is some evidence that the velocity 
dispersion may be somewhat high, but not suffi¬ 
cient to account for a factor of two. An intriguing 
suggestion is that the B gc value is significantly 
low due to the cluster being in an evolutionary 
state different from the others, with more evolved 
galaxies and a very dominant central galaxy. 

5.4. Comparisons with Abell Clusters 

The CNOC1 sample is one of the most homo¬ 
geneous samples for investigating the correlations 
of various properties. The parent sample of the 
CNOC1 sample is the EMSS (with one additional 
Abell cluster with similar X-ray properties), which 
at the time when CNOC1 was carried out was the 
largest and most homogeneous X-ray cluster cat¬ 
alog. Both the optical and X-ray data (Yee et 
al. 1996, and Lewis et al. 1999) were obtained in 
a systematic manner and predominantly from a 
single source. Nevertheless, the CNOC1 sample is 
very small, and it would be instructive to see if a 
larger sample of clusters with data coming from 
more heterogeneous sources follow these correla¬ 
tions. Such a comparison will give confidence on 
the use of these relationship as predictors of vari¬ 
ous cluster properties. 

An excellent set of clusters for such a compari¬ 
son is the LOCOS sample of Lopez-Cruz and Yee 
(see Lopez-Cruz 1997 and YL99). The LOCOS 
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is an imaging survey of 46 low-redshift Abell clus¬ 
ters chosen from the compilation of X-ray detected 
Abell clusters from Jones & Forman (1999). The 
photometric catalogs were created using a proce¬ 
dure identical to that of the CNOC1 survey, with 
the exception that the R c filter band was used in¬ 
stead of Gunn r. We remove the small number 
of Abell richness class 0 clusters from the LOCOS 
samples, as the survey sampled fairly only clusters 
with richness class > 1 (see YL99). Data for the 
velocity dispersion, Tx, and Lx are obtained from 
the literature. 



Fig. 11. - Velocity dispersion versus B gc for Abell 
clusters (open circles). Also plotted for compari¬ 
son are the CNOC1 clusters (solid circles without 
error bars) and the power-law fit from the CNOC1 
clusters. 

For the correlation with velocity dispersion, we 
use data from Fadda et al. (1996) and Girardi et 
al. (1993). YL99 showed a g \ vs B gc plot using all 
clusters with determined using 10 or more ve¬ 
locities, producing a correlation with a relatively 
large scatter. Here we limit our a\ data to those 
comparable to the CNOC1 survey, resulting in 16 
clusters with measured using 25 or more veloc¬ 
ities. The result is shown in Figure 11. Plotted 
for comparison are both the CNOC1 points and 
the best CNOC1 fit. It can be seen that the Abell 


data points fall on the relation on average, albeit 
with a somewhat larger scatter than the CNOC1 
clusters, with the most discrepant clusters being 
at lower richness. The deviation of the Abell clus¬ 
ter oi’s from that predicted by their B gc has a rms 
scatter of 170 km s -1 , and an offset of-55 km s -1 . 
The scatter is equivalent to 17% of the average 
velocity dispersion of the sample. The scatter is 
larger than the CNOC1 sample with MS1455+22 
removed, but slightly smaller than the full CNOC1 
sample. Thus, using a completely different sam¬ 
ple, selected with different criteria and with ve¬ 
locity dispersions obtained from different sources, 
we have found that the B gc values can predict the 
velocity dispersion with an uncertainty of better 
than ~ 20% (about 150 to 200 km s -1 ), and sim¬ 
ilar to that indicated by the CNOC1 sample it¬ 
self. This demonstrates that the richness param¬ 
eter B gc is a useful predictor of the velocity dis¬ 
persion of a cluster, and hence an estimator of its 
mass. The two low-velocity dispersion Abell clus¬ 
ters which show the largest deviations from the 
fitted relation may be an indication of the nature 
of the Abell cluster sample which favors finding 
clusters with projection contaminations. 

For the correlation between richness and X-ray 
properties, we use the compilations of Tx and 
Lx from Wu et al. (1999). Figure 12 plots the 
Tx vs B gc relationship, with the CNOC1 data 
points and best fit indicated. The Abell clus¬ 
ter points fall slightly below the CNOC1 relation, 
with a similar slope and a larger scatter, espe¬ 
cially towards the low Tx direction. However, we 
note that the most discrepant points are also the 
ones with the largest Tx error bars, and most clus¬ 
ters are within 2cr of the regression line from the 
CNOC1 data. 

We find that the Lx vs richness relationship for 
the LOCOS sample has the largest scatter of all 
the correlations discussed. When all LOCOS clus¬ 
ters with available Tx from Wu et al. (1999) are 
plotted, there is only a very rough suggestion of 
a correlation (Figure 13). A closer examination 
shows clearly that the magnitude of the scatter 
in the data is a strong function of the redshift of 
the cluster, with the largest scatter belonging to 
clusters at the low end of the redshift range of 
the sample. On Figure 13, we plot selectively the 
LOCOS data for clusters with z > 0.065. The 
scatter is large, as indicated above, although the 
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Fig. 12.— Tx versus B gc for Abell clusters (open 
circles). Also plotted for comparison are the 
CNOC1 clusters (solid circles without error bars) 
and the power-law fit from the CNOC1 clusters. 

points basically fall along the relation found by 
the CNOC1 sample. The z < 0.065 clusters are 
plotted as crosses, and it is clear that most of 
these objects have very significantly lower Lx than 
that predicted by the CNOC1 relation. The much 
poorer correlation likely indicates the inhomoge¬ 
neous nature of the Lx data compiled by Wu et 
al. (1999). There may also be systematic effects 
between the methods of Lx determination used 
for the CNOC1 sample and the Abell clusters. 
Determining a consistent set of Lx from a het¬ 
erogeneous database is complicated, as data from 
different bands with different resolutions, aperture 
sizes, and signal-to-noise ratios are used. A num¬ 
ber of z < 0.065 clusters have Lx 's which are a fac¬ 
tor of several lower than the expected value, which 
may arise from incorrect aperture corrections due 
to the large angular sizes of these very low-redshift 
clusters. We note that well-studied clusters at low 
redshift, such as Abell 1656 (Coma), in fact do fall 
along the CNOC1 relation. 

For both the Tx and Lx correlations with 
B gc there is a tendency for the Abell sample to 
fall below the CNOC1 regression line, indicating 
a larger B gc relative to Tx or Lx- This difference 



Fig. 13.— X-ray luminosity versus B gc for Abell 
clusters from YL99 with z > 0.065 (open cir¬ 
cles). Also plotted for comparison are the CNOC1 
clusters (solid circles without error bars) and the 
power-law fit from the CNOC1 clusters. Abell 
clusters with x < 0.065 are also plotted, without 
error bars, as crosses. 

can arise from several sources, including simply 
the systematics in the data sets used. One inter¬ 
esting possibility is that this is an indication of the 
relative biases between cluster samples selected us¬ 
ing X-ray and optical criteria. The Abell sample is 
chosen via optical galaxy over-density, and hence 
favors high B gc clusters, while the EMSS sample, 
being selected via X-ray fluxes, preferentially finds 
high X-ray luminosity, hot clusters. Large homo¬ 
geneous samples at similar redshifts and selected 
using both methods are needed to investigate the 
existence of such biases. 

6. Summary 

Using a small but homogeneous sample of clus¬ 
ters from the CNOC1 survey, we have demon¬ 
strated that the optical richness of the clusters, 
measured by the galaxy-cluster center correlation 
amplitude B gc , correlates well with various global 
properties of the clusters. In particular, the corre¬ 
lation of B gc with velocity dispersion is excellent, 
with a rms scatter of only 100 km s -1 when a sin- 
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gle outlier is removed, and 175 km s _1 when the 
whole sample is used. The value of B gc as an esti¬ 
mator of velocity dispersion is verified by using a 
sample of low-redshift Abell clusters, producing a 
similar scatter around the best fit. 

By extension, B gc is correlated with the dynam¬ 
ical mass (as represented by M 200 ) and the virial 
size (as expressed by ?’ 2 oo) of the clusters which 
are derived based on the velocity dispersion. Sim¬ 
ilarly, B gc shows a strong correlation with Tx, and 
a somewhat more scattered correlation with Lx- 
The forms of the correlations of the various prop¬ 
erties with B gc are consistent with scaling laws 
expected based on a simple spherical density dis¬ 
tribution, providing confidence in the existence of 
the underlying correlations. We also compared the 
relative merits of B gc , Tx, and Lx as estimators 
of the virial mass of the clusters. We find that all 
three produce reasonably tight correlations, with 
a scatter of ~ 30%. However, the good correlation 
of Lx with mass is predicated on proper correc¬ 
tions to the Lx determination, which require the 
knowledge of Tx and the spatial distribution of 
flux, making Lx a less useful predictor of mass. 
We note that the correlation of Tx with dynam¬ 
ical mass is consistent with that found by other 
samples using masses determined via X-ray data. 

The correlations investigated here demonstrate 
that we should be able to use optical richness, a 
relatively inexpensive measurement, to estimate 
important attributes of a cluster. This method 
is especially powerful in view of the current capa¬ 
bility of conducting large optical cluster surveys 
with improved cluster finding techniques. It al¬ 
lows samples of tens of thousands of clusters, for 
example, in a 1000 square degree survey, to be 
characterized. An important example of applica¬ 
tions of these scaling relations is the determina¬ 
tion of the mass function of galaxy clusters us¬ 
ing large samples at different redshifts based on 
optical photometric data, being carried out by 
the Red-Sequence Cluster Survey out to 2 ~ 1 
(Gladders 2001; Gladders & Yee 2003). A sim¬ 
ilar method of using total optical light to esti¬ 
mate cluster mass to derive the mass function of 
galaxy clusters at zero redshift has been applied 
by Bahcall et al. (2002) using imaging data from 
the SDSS survey. 

The correlations discussed here are based on 
a small sample of clusters over a relatively wide 


range of redshift of 0.17 < 2 < 0.55. To be able to 
use these scaling relations reliably for clusters over 
a wide range of redshifts, it is important to further 
calibrate these correlations covering a number of 
redshift bins using comparable or larger samples 
from cluster surveys with well-understood selec¬ 
tion functions. These calibrations can be obtained 
in the near future as cluster catalogs from large, 
deep cluster surveys are becoming available. 
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